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ABSTRACT 



Aims. The importance of disc-disc major mergers in galaxy evolution remains uncertain. We study the major merger fraction in a 
SPlTZER/IRAC-selected catalogue in the GOODS-S field up to z ~ 1 for luminosity- and mass-limited samples. 
Methods. We select disc-disc merger remnants on the basis of morphological asymmetries/distortions, and address three main sources 
of systematic errors: (i) we explicitly apply morphological K-corrections, (ii) we measure asymmetries in galaxies artificially red- 
shifted to Zd = 1-0 to deal with loss of morphological information with redshift, and (iii) we take into account the observational errors 
in z and A, which tend to overestimate the merger fraction, though use of maximum likelihood techniques. 

Results. We obtain morphological merger fractions (/m'''') below 0.06 up to z ~ 1. Parameterizing the merger fraction evolution with 
redshift as f^^iz) = /m''''(0)(l + z)'", we find that m = 1.8 ± 0.5 for Mb < -20 galaxies, while m = 5.4 ± 0.4 for > 10'° 
galaxies. When we translate our merger fractions to merger rates (^R™'"''), their evolution, parameterized as = 11l"J''''(0)(l +z)", 

is quite similar in both cases: n = 3.3 ± 0.8 for Mg < -20 galaxies, and n = 3.5 ± 0.4 for Mt, > 10'" Mq galaxies. 
Conclusions. Our results imply that only ~ 8% of today's M-,, > 10'" Mg galaxies have undergone a disc-disc major merger since 
z ~ 1. In addition, ~ 21% of M-^ > 10'" Mq galaxies at z ~ 1 have undergone one of these mergers since z ~ 1.5. This suggests that 
disc-disc major mergers are not the dominant process in the evolution of My, > 10'" Mq galaxies since z ~ 1, with only 0.2 disc-disc 
major mergers per galaxy, but may be an important process at z > 1, with ~ 1 merger per galaxy at 1 < z < 3. 

Key words. Galaxies:evolution — Galaxies:formation — Galaxies:interactions 



1. INTRODUCTION 

The colour-magnitude diagram of local galaxies shows two 
' distinct populations: the red sequence, consisting primarily of 
old, spheroid-dominated, quiescent galaxies, and the blue cloud, 
formed primarilj^y spiral and irregular star-forming galaxies 
.'^ \ jStrateva et all [20011: iBaldrv et al.l 12004). Thi s bimodality has 
, been traced at increasingly higher redshifts (Bell et aU |2004 
up to z ~ 1 : lArn outs et al. 2007; Cirasuolo et al. 200'3, up to 
5^ . Z ~ 1.5; ' Giallong oetalji2005c iCassata et al.l l2008l up to z ~ 2; 
iKriek et al. 2008, at z ~ 2.3). More massive galaxies were the 
first to popula te the red seq uence as a result of the so-called 
"downsizing" jCowie et al.ll996) : massive galaxies experienced 
most of their star formation at early times and are passive by 
z ~ 1, while many of the less massive galaxies have extended star 
formation histories (see Bundv et a l. 2006; S carlata et al..,2007; 
iPerez-Gonzalez et ani2008l and references therein). 

These results pose a challenge to the popular hierarchical 
A-CDM models, in which one expects that the more massive 
dark matter halos are the final stage of successive minor halo 
mergers. However, the treatment of the baryonic component is 
still unclear. The latest models, which include radiative cooling, 
star formation, and AGN and supernova fee dback, seem to re- 
produce the observational trend s better (s ee Bower et al. 2006'; 
Pe Lucia & Blaizot 2007; Stew art et all 2009; Hopkins et al. 
2009bl and references therein). Within this framework, the role 



relative importance in the evolution of galaxy properties, i.e. 
colour, mass, or morphology, is an important open question. 

The merger fraction, /m, defined as the ratio between the 
number of merger events in a sample and the total number of 
sources in the same sample, is a useful observational quantity 
for answering that question. Many studies have determined the 
merger fraction and its evolution with redshift, usually parame- 
terized as /m(z) = /m(0)(l -Hz)'", using different s ample selections 
and methods , such as morpholog i cal cr it eria (IConselice et al.l 
200l l2008l l200St iLaverv et alJ l2004t I Cassata et al.' 12001 



Lotz et alT 



2008i 



Bridge et alj l2007t iKampczvk et al. 200] 



Jogee et al. 2009), kinematic close c ompanio ns (Pattonetajj 
2000 , 2002; Patton & Atfield 2008; iLin et al. 2004, 200i 
De Propris et akf |2005[ l2007t Bluck et al. 2009), sp atiall 



clos e pairs dLeFe vre etal.' '2000; Bundy et al. 2()P, HOO 



1 



of galaxy mergers in the build-up of the red sequence and their 



Brid ge et al.ll2007t Kartaltepe et al. 2007; Hsieh et al. 20081), or 
the correlation function (jBelletal. 2006b; Masjedi et al. 2006!). 
In these studies the value of the merger index m at redshift z < 1 
varies in the range m = 0-4. A-CDM models predict m ~ 2-3 
(Kolatt etal."1999^, 'Go vernatoet aPl 19991: lGottloberetal]l2001l: 
Fakhouri & Ma 2008) for dark matter halos, while suggesting 
a weaker evolution, m ~ 0-2, f or the galaxy merger fraction 
dBerrier et alJl2006l;IStewart et alJi2008) . 

To constrain the role of disc-disc major mergers in galaxy 
evolution, in this paper we study their redshift evolution up to z ~ 
1 in a SPITZER/IRAC-selected catalogue of the GOODS-S area. 
We use morphological criteria, based on the fact that, just after 
a merger is complete, the galaxy image shows strong geomet- 
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rical distortions, particularly asymmetric distortions dConselicd 
1200 3V Hence, high values in the automatic asymmetry index 
A (Abraham et al. 1996; Conselice et al. 2000) are assumed to 
identify disc-disc major merger systems. This methodology 
presents s everal systematic effects, such as signal -to-noise de- 
pendence dConselicd 120031; IConselice et al.l |2005|) or contam- 
ination by non-interacting galaxies with high asymmetry val- 
ues (Jogee et al. 2009; Miller et al. 2008), which lead to biased 
merger fractions if not treated carefully . In a previous study 
of the Groth field, iLopez-Saniuan et alj (l2009l L09 hereafter) 
demonstrated a robust procedure to determine morphological 
merger fractions (/j"''*') using galaxy asymmetries. In that study 
they avoid the loss of information with redshift by artificially 
moving all sources to a common redshift, while the experimen- 
tal error bias, which tends to overestimate the merger fraction 
up to 50%, was addre ssed through use of a maximum likelihood 
method developed in lLopez-Saniuan et aTl (l2008l LGB08 here- 
after). L09 find that the merger rate decreases with stellar mass 
at z = 0.6, and that 20-35% of present-day Mb < -20 galaxies 
have undergone a disc-disc major merger since z ~ 1 . 

This paper is organized as follows: in Sect. |2]we summarize 
the GOODS-S data set that we use in our study, and in Sect.[3]we 
develop the asymmetry index calculations and study their vari- 
ation with redshift. Then, in Sect. |4] we use the methodology to 
obtain the morphological merger fraction by taking into account 
the observational errors. In Sect. |5] we summarize the obtained 
merger fractions and their evolution with z, while in Sect. |6] 
we compare our results with other authors. Finally, in Sect. [T] 
we present our conclusions. We use Hq = 70 kms"' Mpc"', 
£1m = 0.3, and Qa - 0.7 throughout. All magnitudes are Vega 
unless noted otherwise. 

2. DATA 

2.1. The GOODS-S SPITZER/IRAC-selected catalogue 

This work is based on the analysis of the structural param- 
eters of the galaxies catalogued i n the GOODS-Sout h field 
by the Spitzer Legacy Team (see iGiavalisco et al. l 2004). We 
used the Version 1.0 catalogue^ and reduced mosaics in the 
F435W (B435), F6Q6W (Vm), F115W (/ws), and F85QLP (z&so) 
HST/ACS bands. These catalogues were cross-correlated us- 
ing a 1.5" search radius with the GOODS-S IRAC-selected 
sample in the Rainbow Co smologi cal Database p ublished in 
iPerez-Gonzalez et al.l (l2008i see also , Perez-Gonzalez et al..,2005l 
and Barro et al. 2009, in prep.), which provided us with spec- 
tral energy distributions (SEDs) in the UV-to-MIR range, well- 
calibrated and reliable photometric redshifts, stellar masses, star 
formation rates and rest-frame absolute magnitudes. 

We refer the reader to the above-mentioned papers for a more 
detailed description of the data included in the SEDs and the 
analysis procedure. Here, we summarize briefly the main char- 
acteristics of the data set. We measured consistent aperture pho- 
tometry in several U V, optical, NIR and MIR ba nds with the 
method described in iPerez-Gonzalez et all d2008l) . UV-to-MIR 
SEDs were built for aU IRAC sources in the GOODS-S region 
down to a 75% completeness magnitude [3.6]=23.5 mag (AB). 
These SEDs were fitted to stellar population and dust emission 
models to obtain an estimate of the photometric redshift (Zphot), 
the stellar mass (M*), and the rest-frame absolute B-band mag- 
nitude (Mb). 



' http://archive.stsci.edu/prepds/goods/ 



The median accuracy of the photometric redshifts at z < 
1.5 is Izspec - Zphotl/(l + Zspfir) = 0-04, w ith a fraction <5% 
of catastrophic outliers dPerez-Gonzalez et al. 2008, Fig. B2). 
Rest-frame absolute B-band magnitudes were estimated for each 
source by convolving the templates fitting the SED with the 
transmission curve of a typical Bessel B filter, taking into ac- 
count the redshift of each source. This procedure provided us 
with accurately interpolated B-band magnitudes including a ro- 
bustly estimated A:-correction. Stellar masses were estimated 
using exponent ial star formation PEGASEOl models with a 
ISalpeteii (Il955h IMF and various ages, metallicities and dust 
contents included. The typical uncertainties in the stellar masses 
are a factor of ~2 (typical of any stelto popu lation study; see, 
e.g.. iPapovich et al.il2006ilFontana et al.ll2006h . 

Finally, our methodology requires the eiTors in Zphot to be 
Gaussian (Sect. [H LGB08, L09), while Zphot confidence inter- 
vals given by methods do not correlate wi th the differences 
between Zspec's and Zphot's dOvaizu et al.ll2008h . Because of this, 
and following L09, we use 0-,^^^^ = osXl + Zphot) as the Zphot er- 
ror, where cr^_ is the standard deviation in the distribution of the 
variable 6^ = (Zphot - Zspec)/(1 +Zphot), which is well described 
by a Gaussian with mean fig^ ~ and standard deviation cr^„ . We 
found that crg^ increases with redshift, and we took cr^_ = 0.043 
for z < 0.9 sources and o-s_ = 0.05 for z > 0.9 sources. This pro- 
cedure assigns the same error to sources with equal Zphot, but it 
is statistically representative of our sample and ensures the best 
Gaussian approximation of Zphot errors in the merger fraction de- 
termination (Sect. nil. 



2.2. Luminosity- and mass-seiected samples 

The aim of this study is to determine the galaxy merger fraction 
in B-band luminosity- and stellar mass-selected samples. The B- 
band study is motivated by previous studies, which usually se- 
lected their samples in that band. This permits us to compare 
our results with other authors (Sect. 16.2b . Moreover, the stel- 
lar mas s is a fundamental galaxy property that corr elates with 
colour dBaldry et alJl2004l) and morphology dConseli ce 2006a). 

To determine the luminosity limit in the B-band we cal- 
culated the third quartile of the Mb source distribution at 
different redshifts, taking this as a limiting magnitude (e.g., 
iPerez-Gonzalez et all I2OO8I). In the upper panel of Fig. [T] we 
show Mb vs redshift up to Zmax = 1.3 (grey dots) and the lim- 
iting magnitude at different redshifts (black bullets). The upper 
redshift limit in our study, Zmax = 1-3, is fixed by the reliabil- 
ity of the asymmetry index as a morphological indicator with- 
out performing morphological K-corrections (see Sect. 13.1.21 
for details). The black solid curve is the least-squares fit of the 
limiting magnitudes by a third-degree polynomial. At redshift 

Zmax - 1-3, MB^iim 19.5, SO wc sclectcd for our study sources 

with Mb < -19.5. 

We took as limiting mass at each redshift the stellar mass for 
which the IRAC catal ogue i s 75% complete for passively evolv- 
ing galaxies (see lPerez-Gonzalez et alJ2008h . In the lower panel 
of Fig. [T] we show log(M*/Mo) vs redshift up to Zmax = 1-3 
(grey dots) and the 75% of completeness at different redshifts 
(black bullets). The black solid curve is the least-squares fit of 
the completeness points by a power-law function. At redshift 
Zmax - 1-3, log(M*jin,/MQ) ~ 9.8, SO we selected sources with 

> 10'" Mq for our study. 
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Fig. 1. Top: distribution of Mb vs redshift for IRAC catalogue 
sources. The black dots are the limiting magnitude of the sur- 
vey at each redshift, defined as the third quartile in magni- 
tude distributions. The solid black curve is the best fit of the 
limiting magnitude points by a third-degree polynomial. The 
black dashed line shows the Mb - -19.5 limit of our study. 
Bottom: distribution of log(M*/MQ) vs redshift for IRAC cat- 
alogue sources. The black solid curve shows the stellar mass 
above which the sample is 75% comple te for passively evolv- 
ing galaxies jPerez-Gonzalez et aT Mm . The black dashed line 
shows the log(M*/M0) = 10 limit of our study. 

3. ASYMMETRY INDEX 

The automatic asymmetry index (A) is one of the CAS morpho- 
logical indices (IConselicell2003h . This index is defined as 



2 I/O - /i8ol XlBo-B 



1801 



Z |/ol 



2 |/ol 



(1) 



where Iq and Bq are the original galaxy and background images, 
/i80 and Bi8() are the original galaxy and background images ro- 
tated 180 degrees, and the summation spans all the pixels of the 
images. The background image is defined in detail in the next 
section. For further d etails on the asymmetry calculation see 
IConsehce et alJ (l2000l) . This index gives us information over the 
source distortions and we can use it to identify recent merger sys- 
tems that are highly distorted. In previous studies a galaxy was 
taken to be a recent m erger if its asymmetry index is A > A^ , 
with Am = 0.35 ( e.g., IConselicd 120031: iDe Propris et"ani2007t 
[Bridge et alj|2007|) . This methodology presents several system- 
atic effects, such as sig nal-to-noise dependence (Conselice 2003t 
IConselice et alj|2b05l) . contamin ation by non-interacting galax- 
ies with high asymmetry values jJogee et alj|2009l : iMiller et"an 



2008|), contamination by nearby bright sources dPe Propris et al.l 

2007) , or the pass-band in which we measure the asymmetry 
(Cassata et al. 2005; Tavlor-Mager et al. 2007; Conselice et aT| 

2008) , which must be carefully treated to avoid biased merger 
fractions. In the following sections we detail how we determined 
the asymmetry index and its dependence on several factors, such 
as the background image Bq that we use (Sect. 13.1. It , the pass- 
band in which we calculate it (Sects. [Tl. 2113. 1.3T l and the signal- 
to-noise of the source (Sects. [3?1.4ll3.2|[33] l. 

3.1. Asymmetry calculation 

3.1.1. Background dependence 

In Eq. ([U) we have a dependence on the background image Bo; 
that is, different backg round images yield di fferent asymmetries 
for the same source ( Conselice et alj|2003h . To minimize this 
effect we determined the asymmetry of each source with five 
different background images. These background images are sky 
source-free sections of 50x50 pixels located in the same posi- 
tion in the four HST/ACS filter images, and were chosen to span 
all the GOODS-S area. The asymmetry of one source was the 
median of those five background-dependent asymmetries. 

3.1.2. Pass-bands and redshift range 

Galaxy morphology depends on the band of observation (e.g. 
Kuchinski etal] 120001; iLauger et al.1 120051; iTavlor-Mager et al.l 
20071) . In particular, when galaxies contain both old and young 
populations, morphologies may change very significantly on 
both sides of the Balmer/4000A break. The asymmetry index 
limit Am - 0.3 5 was established in the rest-frame B-band 
dConselic^ l2003h . When dealing with galaxies over a range of 
redshifts, in order to avoid systematic passband biases with red- 
shift, one needs to apply a so-called morphological K-correction 
by performing the asymmetry me asurements in a ban d as close 
as possible to rest-frame B (e.g.. ICassata et alll2005b . or apply 
statistical corrections for obtaining asymmetries in rest-frame B 
from asymmetry measurements in rest-frame U dConsehce et al.l 
2008). Taking advantage of the homogeneous multiband imag- 
ing provided by the GOODS survey, we entirely avoid morpho- 
logical K-correction problems in the present study by perform- 
ing asymmetry measurements on all GOODS-S B435, Veoe, '775, 
and Z850 images, and using for each source the filter that most 
closely samples rest-frame B. 

To determine the redshift ranges over which rest B-band or 
f/-band dominates the flux in the four observational HST/ACS 
filters, B435, V(,o6, '775, and zs,5q, we defined the function 



fRF(z) = 



(2) 



where Prf and Pacs are the transmission curves of the rest- 
frame reference filter and one HST/ACS filter, respectively. In 
Fig.|2]we show the function /b(z) for the four ACS filters (black 
curves), and fu(z) for zgso (grey curve). On the basis of this fig- 
ure, B435 asymmetries were used for < z < 0.15 sources; Vgoe 
asymmetries for 0.15 < z < 0.55; ijjs for 0.55 < z < 0.9; and 
Z850 for 0.9 < z < 1.3. Staying within rest-frame B imposed a 
maximum redshift of Zmax = 1.3. 

Note that, because the ML method used in the merger frac- 
tion determination (Sect. |4]i takes into account the experimental 
errors, we had to include in the samples not only the sources 
with Zi < Zup, where Zup is the upper redshift in our study, but 
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Fig. 2. Function for the four ACS filters: B435 (black dashed 
curve), Vdob (black dotted curve), i-n^ (black dot-dashed curve), 
and Z850 (black solid curve). The grey solid curve is the func- 
tion fuiz) for the z^so filter. The vertical black solid line is the 
maximum redshift, Zmax - 1-3, in our study. 

also sources with z, - 2cri < Zup in order to ensure completeness. 
Because of this, Zup must fulfil the condition Zmax - 2cr^,(l -H 
Zmax) = Zup, which yields Zup ~ 1.1- We took as minimum red- 
shift in our study Zmin -0.1 because of the lack of sources at 
lower redshifts. This yields Zdown = Zmin + 20-^^(1 + Zmin) ~ 0.2, 
which ensures completeness and good statistics. Applying these 
redshift limits we finally have 1740 galaxies with Mb < -19.5 
and 982 with > lO'*' Mq. The number of galaxies quoted 
here was obtained after removing problematic border sources 
(Sect, mi- 

3.1 .3. Determining tine asymmetry of sources witli 
pliotometric redsliifts 

Roughly ~ 40% of the sources in our samples do not have 
spectroscopic redshifts and we rely on photometric redshift de- 
terminations. In these cases, our source could have its rest- 
frame B-band flux in two observational ACS filters, within Icr. 
To take this into account we assumed three different redshifts 
for each photometric source: z^^^^^ = Zphot - cr;^,,,^,, Zphot, and 
■^phot ~ -^phot + Czphot ■ We determined the asymmetry in these three 
redshifts. We then performed a weighted average of the three 
asymmetry values such that: 

Ao = 0.16A(Zpi,„,) + 0.16A(z;,„.) + 0.68A(zphot), 

where A(z) is the asymmetry of the source at redshift z- We used 
the same average procedure with the uncertainties of the three 
asymmetries and added the result in quadrature to the rms of the 
three asymmetry values to obtain cr^„ . In sources with Zspec we 
only determined the asymmetry at the source redshift. L09 show 
that the two different asymmetry determinations do not intro- 
duce systematic differences between sources with and without 
spectroscopic information. 

3.1 .4. Boundary effects and briglit source contamination 

The signal-to-noise in HST/ACS decreases near the boundaries 
of the images, where the exposure time is lower. This affects 
our asymmetry values in two ways: the SExtractor segmenta- 
tion maps that we use to calculate the asymmetry have many 
spurious detections, and any of the five backgrounds defined in 
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Fig. 3. Asymmetry vs redshift in the Mb < -20 sample (grey 
dots in both panels). Top: asymmetries of the sources measured 
on the original images. Bottom: asymmetries of the sources mea- 
sured on images artificially redshifted to Zd = 0.1. Open squares 
in both panels are the mean asymmetries in 0.1 redshift bins. 
The black solid line is the least-squares linear fit to the mean 
asymmetries in the [0.5,1.3) redshift interval. 



Sect. 13.1.11 is representative of the noisier source background. 
Th e problem with segmen tation maps was noticed previously 
by iDe Propris et al.l (l2007h . where the segmentation maps for 
50% of their initial 129 galaxies with A > 0.35 are incorrect, 
or are contaminated by bright nearby sources. With this in mind, 
we visually inspected all the sources looking for boundary or 
contaminated sources. We found that boundary sources had sys- 
tematically high asymmetry values, and had segmentation maps 
contaminated by spurious detections. To avoid biased merger 
fraction values we excluded all border sources (high and low 
asymmetric) from the samples. We found only two sources con- 
taminated by bright nearby sources. For these we redefined the 
SExtractor parameters to construct correct segmentation maps 
and redetermined the asymmetry. 

3.2. Asymmetries at a Reference Redshift 

The asymmetry index measured on survey images systematically 
varies with the source redshift due, first, to the (1 -Hz)"* cosmolog- 
ical surface brightness dimming, which can modify the galaxy 
area over which asymmetry is measured, and, second, to the loss 
of spatial resolution with z. Several papers have attempted to 
quantify these effects by degrading the image spatial resolution 
and flux to simulate the appearance that a give n galaxy would 
have at different redshifts in a given survey. IConselice et al.l 
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( I2003L I2008h : and ICassata et al] (l2005h degraded a few local 
galaxies to higher redshifts and found that asymmetries decrease 
with z. [Conselice et al. (2003) also noted that this decrease de- 
pends on image de pth, and that luminous galaxies are less af- 
fected. In addition, IConselice et alj (l2005i) show that irregular 
(high asymmetry) galaxies are more affected than ellipticals (low 
asymmetry). A zeroth-order correction for such biases was im- 
plemented by Conselice et al. (2003, 2008, 2009) who applied a 
AA^ term, defined as the difference between the asymmetry of lo- 
cal galaxies measured in the original images and the asymmetry 
of the same galaxies in the images degraded to redshift z- Their 
final, corrected asymmetries are Af = Aq + AAr, where Aq is the 
asymmetry measured in the original images. With these correc- 
tions, all the galaxies have their asymmetry referred to z = 0, 
and the local merger criterion A > A^ = 0.35 is then used. 

In their study, L09 improve on the above procedure, and we 
apply their methodology to our data set. We compute a correc- 
tion term individually for each source in the catalogue, but rather 
than attempting to recover z - values for A we degrade each 
of the galaxy images to redshift Zd = 1 ; we then obtain our final 
asymmetry values Af directly from the degraded images. With 
this procedure, we take into account that each galaxy is affected 
differently by the degradation; e.g. the asymmetry of a low lu- 
minosity irregular galaxy dramatically decreases with redshift, 
while a luminous elliptical is slightly affected. We choose Zd = 1 
as our reference redshift because a source at this (photometric) 
redshift has Zd + o-,^ ~ Zup = 1 . 1 ; that is, the probability that our 
galaxy belongs to the range of interest is ~ 85%. Because we 
work with asymmetries reduced to Zd - 1, the asymmetry crite- 
rion for mergers, A^, needs to be reduced to z = 1 . We discuss 
this in Sect. [33] 

We have already mentioned that ~60% of the sources in 
the samples have spectroscopic redshifts, hence redshift infor- 
mation coming from photometric redshifts for the remaining 
~ 40% of the sources has large uncertainties. As in the Aq cal- 
culation process (Sect. 13.1.31 Eq. Q), to take into account the 
redshift uncertainty when deriving the asymmetries at Zd = 1 
we started from three different initial redshifts for each source, 
Vot " ^ph°t - «"<-pho.' ^phot, and z+ho, = Zphot + o-zp„„„ and degraded 
the image from these three redshifts to Zd - 1. We then per- 
formed a weighted average of the three asymmetry values such 
that 

Af = 0.16Ai(Zphot) + 0.16Ai(z;hot) + 0.68Ai(zphot), (4) 

where Ai(z) denotes the asymmetry measured in the image de- 
graded from z to Zd = 1. When a spectroscopic redshift was 
available, the final asymmetry was simply Af = Ai(zspec)- We 
did not apply any degradation to sources with z > 1 ; that is, we 
assumed that Ai(z > 1) = Aq. 

To obtain the error of the asymmetry, denoted by 0-^^, for 
sources with photometric redshifts, we averaged the uncertain- 
ties of the three asymmetries following Eq. © and added the re- 
sult in quadrature to the rms of the three asymmetry values. The 
first term accounts for the signal-to-noise error in the asymmetry 
value, while the second term is only important when differences 
between the three asymmetry values cannot be explained by the 
signal-to-noise first term. In sources with spectroscopic redshifts 
we took as cr^, the uncertainty of the asymmetry A i(zspec)- 

The degradation of the images was performed with cos- 
MOSHiFT (iBalcelLs et al.ll20()3l) . which performs repixelation, psf 
change and flux decrease over the sky-subtracted source image. 
The last cosmoshift step is the addition of a random Poisson sky 
noise to the degraded source image to mimic the noise level of 



Table 1. Degradation rate for different luminosity samples 



Sample selection 


'^tot 


Sa 




Ao 




(1) 


(2) 




(3) 


Mb < -19.5 


1740 


-0.09 ± 


0.01 


0.156 


Mb < -19.75 


1402 


-0.07 ± 


0.01 


0.160 


Mb < -20 


1122 


-0.05 ± 


0.01 


0.163 


Mb < -20.25 


869 


-0.05 ± 


0.02 


0.165 


Mb < -20.5 


648 


-0.04 ± 


0.02 


0.161 



NOTES. Col. (1) Number of sources in the sample with 0.1 < z < 1.3. 
Col. (2) Degradation rate of the asymmetry, AA = Sa^Z- Col. (3) 
Median asymmetry of sources with z < 1. 

the data. As a result of this last step, two cosmoshift degrada- 
tions of the same source yield different asymmetry values. We 
took the asymmetry of each degraded source, Ai(z), to be the 
median of asymmetry measurements on five independent degra- 
dations of the original source image from z to Zd = 1 . With all 
the aforementioned steps, each Ai(z) determination involved 25 
asymmetry calculations, while the uncertainty in Ai(z) was the 
median of the five individual asymmetry errors. 

The asymmetries Af referred to Zd = 1 provide a homo- 
geneous asymmetry set that permits consistent morphological 
studies in the GOODS-S field (Lopez-Sanjuan et al. 2009b, in 
preparation). 

3.3. Asymmetry Trends with Redshift 

For a sample of galaxies over a range of redshifts, the statistical 
change with z of the measured asymmetries Aq is the combined 
effect of loss of information (as shown in the previous section) 
and changes in the galaxy population. In contrast, the redshift 
evolution of Af reflects changes in the galaxy population alone, 
given that the morphological information in the images used to 
determine Af is homogeneous for the sample. As already dis- 
cussed in L09 for the Groth field, we show here that the z trends 
of Ao and Af are quite different. 

In the top panel of Fig. [3] we show the variation of Aq with 
redshift in a Mb < -20 selected sample, while in the bottom 
panel we see the variation of Af for the same sample. In both 
panels, open squares are the median asymmetries in A(z) = 0.1 
redshift bins, and the black solid line is the best linear least- 
squares fit to the 0.5 < z < 1.3 points. Aq is seen to decrease with 
redshift, A() - 0. 19-0.049z, while theAf distribution is flat,Af ~ 
0.14. For Ao, the negative slope reflects the fact that the loss of 
information with redshift (negative effect on A) dominates over 
genuine population variations (a positive eff ect because galaxies 
at higher redshift are more asymmetric; e.g. lCassata et al.ir2005l; 
IConselice et al.ll2005l) .In Af the information level does not vary 
with the redshift of the source, so we only see population effects. 
In this case the slope is null, but this is a field-to-field effect: L09, 
with the same methodology and sample selection, obtain Af cc 
0.05z. This indicates that we cannot extrapolate results from one 
field to another, and that individual studies of systematics are 
needed. We take as degradation rate (Sa) the difference between 
both slopes and assume that the merger condition A^ varies with 
redshift as An,(z) = Am(0) - 6aZ = 0.35 - 6aZ. 

Is the degradation rate the same for all luminosity selections? 
We expect less asymmetry variation with redshift in bright sam- 
ples, because the y are less affected by cosmological dimming 
(IConselicd l2003l) . We repeated the previous analysis with dif- 
ferent Mb selection cuts, from Mb < -20.5 to Mb < -19.5 
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Fig. 4. Degradation rate 6a vs the selection magnitude Mb of the 
sample (black squares). The black solid line marks 5a - -0.05, 
the estimated degradation rate for Mb < -20 galaxies. 



160 
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?> 140 




Fig. 5. Mb distribution of M* > lO'" Mq galaxies. The black 
solid Une is a Gaussian with ju = -20.5 and cr - 0.9. 



(the latter is the limiting magnitude in our study. Sect. 12.21 ). We 
summarize the results in Table[T]and Fig.lH asymmetry is more 
affected by redshift changes in less luminous samples, as ex- 
pected. Interestingly, the degradation rate is roughly constant up 

to Mb - -20, 6a 0.05 (black soUd Une in Fig. HI, but then 

becomes more pronounced by a factor of 2, 6a 0.09, in only 

0.5 magnitudes. One could argue that the sharp increase of 6a 
for samples including Mb > -20 sources arises because such 
sources have higher initial asymmetry Aq : a faint irregular galaxy 
is more affected by loss of information than a bright elliptical. 
However, we see in the last column of Table [T] that the mean 
asymmetry of sources with z < 1.0 is similar in all samples, 
Ao ~ 0.16. Hence, the degradation rate increases because faint 
sources have lower signal-to-noise than luminous ones. Because 
of this, we decided to restrict our study to the 1 122 sources with 
Mb < -20 to ensure that degradation affects all the galaxies 
in our sample in the same way, making the merger condition 
Am(l) = 0.35 - 6a = 0.30 representative. 

How important is this luminosity dependence for the mass- 
selected sample? The Mb distribution of > 10'" M© galaxies 
is well described by a Gaussian with fi - -20.5 and cr - 0.9, 
Fig. |5] We found that 70% of the galaxies have Mb < -20, and 
that the degradation rate for the whole sample is 6a - -0.05. 
This tells us that the faint sources in this sample do not signifi- 
cantly affect the degradation rate, making the Mb < -20 merger 



condition representative also for the mass-selected sample. In 
conclusion, we used Ain(l) - 0.30 for both samples. 



4. MERGER FRACTION DETERMINATION 



Following IConselic3 (l2006bh . the merger fraction by morpho- 
logical criteria is 



■mph 



"tot + (k- l)n^' 



(5) 



where is the number of the distorted sources with A > Am, 
and «tot is the total number of sources in the sample. If a- > 2 we 
obtain the galaxy merger fraction, f^^, the fraction of galax- 
ies undergoing mergers, and k represents the average number of 
galaxies that merged to produce one distorted remnant. If /c = 1 
we obtain the merger fraction, : the number of merger events 
in the sample. We use k = I throughout this paper. 

The steps we followed to obtain the merger fraction are 
described in detail in LGB08. In this section we provide a 
short summary. If we define a two-dimensional histogram in the 
redshift-asymmetry space and normalize this histogram to unity, 
we obtain a two-dimensional probability distribution defined by 
the probability of having one source in bin [zk,Zii+\) n [A;, A;+i), 
namely p^i, where the index k spans the redshift bins of size Az, 
and the index I spans the asymmetry bins of size AA. We con- 
sider only two asymmetry bins split at A^, such that the prob- 
abilities pti describe highly distorted galaxies (i.e. merger sys- 
tems), while the probabilities p^o describe normal galaxies. With 
those definitions, the morphologically based merger fraction in 
the redshift interval [zk,Zii+\) becomes 



/.mph 



Pkl 



PkO + Pkl 



(6) 



In LGB08 they describe a maximum likelihood (ML) method 
that yields the most probable values of pki taking into account 
not only the z and A values, but also their experimental errors. 
The method is based on the minimization of the joint hkehhood 
function, which in our case is 



L(zi,Ai\p'i^,,o-,.,crA,) 



In ] Z Z T"ERF(z, /, ^)ERF(A, /, 



where 



ERF(77, i, k) = erfj 



V2cr,, 



■ erf 



iq - rik 
V2o-„, 



(7) 



(8) 



In the above equations, &rf{x) is the error function; Zi and A, are 
the redshift and asymmetry values of source /, respectively; cr^, 
and (Ta^ are the observational errors in redshift and asymmetry 
of source /, respectively; and the new variables = \n{pki) 
are chosen to avoid negative probabilities. Equation (|7]l was ob- 
tained by assuming that the real distribution of galaxies in the 
redshift-asymmetry space is described by a two-dimensional 
distribution pki = exp(/9y, and that the experimental errors are 
Gaussian. Note that changing variables to p'^.^ - \n{pki), Eq. ^ 
becomes 



/.mph 

J m,k 



(9) 
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Fig. 6. Asymmetry vs log(M*/MQ) for light-weight (open cir- 
cles) and faint (bullets) samples (see text for details). The black 
solid line shows the merger criterion Am(l) = 0.30. 



LGB08 show, using synthetic catalogues, that the experimen- 
tal errors tend to smooth an initial two-dimensional distribu- 
tion described by pi^i, due to spill-over of sources to neigh- 
bouring bins. This leads to a ~ 10 - 30% overestimate of the 
galaxy merger fra ction in typical observational cases. L09 and 
iLotz et al] (l2008ah find similar trends in their study of the mor- 
phological merger fraction in the Groth Strip. LGB08 addition- 
ally show that, thanks to the use of the ML method, they can 
accurately recover the initial two-dimensional distribution: the 
fractional difference between the input and ML method merger 
fractions is a tiny ~ 1% even when the experimental errors are 
similar to the bin size. That is, the ML results are not biased by 
the spill-over of sources to neighbouring bins. 

We obtained the morphological merger fraction by applying 
Eq. (|9]l using the probabilities recovered by the ML method. 
In addition, the ML method provides an estimate of the 68% 
confidence intervals of the probabilities which we use to 

obtain the f'"f^ 68% confidence interval, denoted [cr"„,„, , cr+ , ] . 

This interval is asymmetric because f^f^ is described by a log- 
normal distribution due to the calculation process (see LGB08 
for details). Note that, in LGB08, k - 2 is used in Eq. 0, but 
the method is valid for any k value. 

We also determined the morphological merger fraction by 
classical counts, f™''!' = nm^^Vntaf^^ where n5^!^^'* is the num- 

' m.class m / tot ' m 

ber of galaxies in a given bin with Af > A, and n'^ll'^^ is the total 
number of sources in the same bin. We obtained the f"''!' un- 

•' m.class 

certainties assuming Poissonian errors in the variables. 

Finally, and following L09, sect. 4.1, we performed simula- 
tions with synthetic catalogues to determine the optimum bin- 
ning in redshift for which the ML method results are reliable. 
The simulations were made in the same way as in L09, so here 
we only report the results of the study: we can define up to 
three redshift bins, namely zi - [0.2, 0.6), zi = [0.6, 0.85), and 
Z3 = [0.85, 1.1). The first bin is wider than the other two, 0.4 vs 
0.25, because of the lower number of sources in the first inter- 
val. In the next section we study the merger fraction evolution 
with redshift with these three bins (§ 15.11 ). We will also provide 
statistics for the zo = [0.2, 1.1) bin in order to compare the ML 
and classical merger fraction determinations. 



5. RESULTS 

We summarize in Table |2] the main characteristics of the two 
samples under study; i.e. the total (ntot) and distorted («„,) num- 
ber of sources, both for classical counts («'^'^'*^) and the ML 
method (n'^^), and major merger fractions. Note that the number 
of ML method galaxies is not an integer. Indeed, the ML method 
gives us a statistical estimate of the probability p^i = exp(/:i^,) 
of finding one source in the redshift bin k, and in the asym- 
metry bin /, so the estimated number of galaxies in that bin, 
W/U.ML - NtotPki^zAA, where A^tot is the total number of galax- 
ies in the sample, need not be an integer. The merger fraction by 
the ML method is roughly half that in the classical determina- 
tion (0.035 vs 0.077 in the luminosity-selected sample, 0.025 vs 
0.050 in the mass-selected sample). This highlights the fact that, 
whenever the spill-over eff'ect of large measurement errors is not 
taken into account, morphological merger fractions can be over- 
estimated by a factor of ~ 2. We use this result later in Sect. 16.31 
and in the next section we use only merger fractions obtained by 
the ML method. 

We find that correction of redshift-dependent biases is 
equally important. If we use the raw asymmetry values deter- 
mined on the original images, and apply the local Universe 
merger selection criterion Ao > 0.35, the resulting merger frac- 
tions come up a factor 2 higher than the ones listed in Table |2] 
Recall that the latter come from Af values homogenised to a 
common reference Zd = 1 (Sect. 13.2b . This emphasises that 
published merger fractions which do not work with redshift- 
homogeneous data, may be significantly biased. Interestingly, 
an identical comparison to the one just described, applied to 
Groth strip data, lead L09 to conclude that redshift effects are 
not important for merger fraction determinations. The different 
behaviour of the Groth data from L09 and our GOODS-S data 
might be due to cosmic variance, or to depth differences between 
the two data sets. In general though, artificial redshifting of the 
galaxies is needed to ensure reliable results. 

Table |2] shows that the merger fraction from the mass- 
selected sample is lower than that from the luminosity-selected 
sample. What is the origin of this difference? To answer this 
question, we define two subsamples: the faint sample (galax- 
ies with Mb > -20 and M* > 10'° Mq), and the light-weight 
sample (sources with M* < lO'" and Mb < -20). The 
faint sample comprises 272 sources, while the light-weight sam- 
ple comprises 408 sources. In Fig. |6]we show both samples in 
the mass-asymmetry plane: light-weight galaxies have higher 
asymmetry, Af = 0.14, while faint galaxies are more symmet- 
ric, Af - 0.07. The light-weight sample comprises 43 sources 
with Af > 0.30 (10.5% of the sample), while the faint sam- 
ple comprises only seven distorted sources (2.5% of the sam- 
ple). These numbers suggest, in agreement with L09, that: (i) 
an important fraction of the B-band high asymmetric sources 
are low-mass disc-disc merger systems that, due to merger- 
triggered star-formation, have their B-band luminosity boosted 
by 1.5 magnitudes (Bekki & Shiova 2001), enough to fulfil our 
selection cut Mb < -20; and (ii) the faint objects are earlier 
types dominated by a spheroidal component which, when sub- 
ject to a major merger, does not distort enough to be picked up 
as merger systems by our asymmetry criterion. 

5.1. Merger fraction evolution 

We summarize in Table [3] the morphological merger fraction at 
different redshifts in GOODS-S. We obtain low merger fractions, 
always lower than 0.06, similar to the L09 results for the Groth 
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Table 2. Sample characteristics in the 0.2 < z < 1.1 range 



Sample selection 


class 
"tot 

(1) 


class 

(2) 


/■mph 
m,class 

(3) 


"tot 

(4) 


(5) 


/■mph 
/m.ML 

(6) 


Mb < -20 


793 


61 


0.077 ± 0.010 


881.9 


30.7 


035+" "'" 


> 10'° Mo 


759 


38 


0.050 ± 0.008 


819.3 


20.2 


025+'' °''^ 



NOTES. Col. (1) Number of galaxies with 0.2 < z < 1.1 by classical counts. Col. (2) Number of distorted galaxies with 0.2 < z < 1.1 and 
Af > 0.30 by classical counts. Col. (3) Morphological major merger fraction by classical counts. Col. (4) Number of galaxies with 0.2 < z < 1.1 
by ML method. Col. (5) Number of distorted galaxies with 0.2 < z < 1.1 and Af > 0.30 by ML method. Col. (6) Morphological major merger 
fraction by ML method. 

Table 3. Morphological major merger fractions /™'''' in GOODS-S 



Sample selection 


z = 0.4 


z = 0.725 


z = 0.975 


f:r(0)-' 


m 


Mb < -20 
M* > 10'° Mo 


023+" "-'^ 
006+°'"'* 


022+"°'' 


043+""'^ 

037+" °"' 
"•"■"-noil 


0.013 ± 0.003 " 
(1.0 ± 0.2) X 10-3 


1.8 ±0.5 " 
5.4 + 0.4 



" Best f7\z) = .C (0)(1 + zY" fit to the dat a. 
This fit includes the lDe Propris et aTl i2007b local value. 

Table 4. Morphological merger fraction in GOODS-S at 0.6 < z < 0.85 



Sample selection 


«LSS 


w/LSS 


w/o LSS 


LSS (z = 0.735) 




(1) 


(2) 


(3) 


(4) 


Mb < -20 


72 


031+""'" 

-QQU 


0.026:";°^° 


044+" "^'^ 
-Q-QIS 


M. > 10'° Mo 


94 


022+" °" 


0.018i»l^ 


032+°-"" 



NOTES. Col. (1) Number of galaxies in the Large Scale Structure (LSS). Col. (2) Merger fraction in the sample with LSS. Col. (3) Merger fraction 
in the sample without LSS. Col. (4) Merger fraction in the LSS. 

field. The merger fraction increases with redshift in both the 
luminosity- and the mass-selected samples, but this growth is 
more prominent in the mass-selected sample. We can parameter- 
ize the merger fraction evolution as 

f7\z) ^ f7\o)(i + zT (10) 

and fit our data. Note that, in the luminosity-selected sample, 
we also use the Mb < -20 estimation from L09 of the Mg < 
-19 local merger fraction, drawn from the MGCQ (Millenium 
Galaxy Catalogue), from lDe Propris et alJ (l2007h : /Z\0.07) = 
0.014+g gy3. Wc summarfzc the results in Tableland Fig.]?] The 
merger index m is higher (3cr) in the mass-selected sample (bul- 
lets) than in the luminosity-selected sample (open triangle for 
[Dc Propris et al. 2007 local value; open squares for our data), 
5.4 vs 1.8, while the merger fraction in the local universe is 
lower in the mass-selected sample, 0.001 vs 0.013. The fact that 
the high er m, the lower f^'^^jO), was predicted by semianalytical 
models toochfar & Burkertll200l . We compare these values 
with those from previous studies in Sect. 16.21 

5.2. Large Scale Structure effect 

It is well known that the more prominent large scale structure 
LSS) in the GOODS- S field is located at redshift z = 0.735 
Ravikumar et al.l l2007!). In order to check the effect of this LSS 



2 www.eso.org/ jliske/mgc 



on our derived merger fractions, we recalculated them by ex- 
cluding the s ources wit hin dv < 1500 km s'' (dz ~ 0.01) of 
z = 0.735 (Rawa tet aDl2008h . In Table H we summarize the 
number of sources in the LSS for each sample (hlss), and the 
previous and recalculated merger fractions, both in the field and 
in the structure. The merger fraction is higher in the LSS than 
in the field. Note that the variation in the field values is well 
reported by the error bars. How does this LSS affect the previ- 
ously inferred merger evolution? If we again fit the data without 
LSS, we find that /m'''^(0) does not change, while the value of 
m decreases only by 0. 1 in both the luminosity- and the mass- 
selected samples, so our conclusions remain the same. We shall 
therefore use the fit values in Table |3]in the remainder of the pa- 
per We concentrate on the LSS at z = 0.735, and ignore other 
structures in GOODS-S. The next two more important ones are 
located at z = 0.66 and z = 1.1. The former is an overdensity 
in redshift space, but not in the sky plane, while the latter is a 
cluster, but comprises an order of ma gnitude fewer so urces than 
the z = 0.735 structure (145 vs 12. iAdami et al.ll2005h . 



6. DISCUSSION 

First we compare our results with merger fraction determina- 
tions from other authors. In Fig. [T] we show our results (open 
squares for Mb < -20 galaxies and bullets for > 10'° Mq 
galaxies). The other points are those from the l i teratu re: the 
Mb < -20 estimate by L09 of the lDe Propris et al.l (l2007b Mb < 
-19 merger fraction; the merger fraction for Z?-band luminosity 
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selected galaxies in AEGIl^ (All- Wavelength Extend ed Groth 
Strip International Survey) from iLotz et aEI d2008 a'): the re- 
sults from Conselice et al. (2009) in COSMO^I (Cosmological 
Evolution Survey) and AEGIS for > 10'° Mq galax- 
ies; and the merger fraction for > 5 x 10'" M© galax- 
ies in GEMlj^ (G a laxy Evolution from Morphology and SEDs) 
from Jogee et al. (2009). Note that the mass selection from 
Jogee et al. (2009) has been adapted to a Salpeter IMF (Salpeter 
19551). All the pr evious merger fractions except those from 
Jogee et a l.' ('2009') are from automatic indices for major merg- 
ers. The Jogee et al. (2009) results are by visual morphology 
and reflect major+minor mergers; the dashed rectangle marks 
their expected major merger fraction. For luminosity-selected 
samples (open sy mbols) our values are in good agre ement with 
iDe Prop ris et al.l ( l2007i) . but are lower than those from lLotz et alj 
(2008a), who apply diff'erent sample selection and merger crite- 
ria from ours and do not correct the effect of observational errors, 
thus making comparison difficult. 

In the mass-selected case our results are in goo d agreement 
with th e expected visual major merger fraction from lJogee et alj 
( |2009 |) (dashed lines), supporting the robustness of our method- 
ology for obtaining major merger fracti ons statistically. Our val - 
ues are significantly lower that those of IConselice et al.l (l2009h . 
especially at z > 0.7, where there is a factor 3 difference . 
The asymmetry calculation performed bv' Conselice et al.l (l2009l) 
does not take into account the spill-over effect of observational 
errors in their merger fraction determination. We show here that 
such effects ma y lead to the higher value obtained by them. 
IConselice et al.l 32009) assume two main statistical corrections 
at z > 0.7: the information degradatio n bias (AA^. § | 3.2[l and the 
morphological K-correction (AA^, see IConselice et al.ll2008l for 
details). The first correctio n is AA^ = 0.5 and has an associated 
uncertainty of ctaa^ ~ 0.08 (IConselice et al.ll2003L Table 1). The 
morphological K-correction depends on redshift; to simplify the 
argument, we do not consider its uncertainty in the following. In 
addition, each source asymmetry has its own signal-to-noise un- 
certainty, which in our study is ~ 0.03 at the se redshifts. We 
reproduced the same methodology applied by IConselice et alj 
(2009) on synthetic catalogs created as in Sect. |4] For further 
details about simulation parameters and assumptions, see L09. 
In the simulations we defined two redshift intervals, namely 
Z2 = [0.6, 0.85) and Z3 - [0.85, 1.1), taking our results in these 
redshift intervals as input merger fractions, f^'^^ = 0.022 in the 
first interval, and - 0.037 in the second. We then extracted 
2000 random sources in the redshift-asymmetry plane, applying 
an asymmetry error to them of cta = ( 3.08, which is represen - 
tative of the asymmetry uncertainties in IConselice et al.l (l2009l) . 
We assumed cr- - for simplicity. Merger fractions wer e de- 
rived from classical histograms as in ' Conselice et al.l(l2009l) . We 
repeated this process 100 times and averaged the results. This 
process yields f^''^ ~ 0.11 in the first interval, an d /^''^ ~ 0.13 
in the second, which is similar to lConsehce et al.l 12009) results 
at these redshifts. In contrast, the ML method was able to re- 
cover the input merger fractions. The exercise demonstrates that 
the observed differences betwen the t wo studies can b e natu- 
rally explained as a bias introduced in IConselice et al.l (l2009l) 
by not accounting fo r spill-over of sources due to observational 
errors. The fact that IConselice et all (l2009l) study is performed 
over ~ 20000 galaxies, 20 times more sources than in our study. 



^ http://aegis.ucolick.org/ 
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cannot correct the errors. As emphasized by LGB08, experimen- 
tal systematic errors are not cured by increasing sample size: the 
ML method is needed. 

6. 1 . Groth vs GOODS-S merger fractions: cosmic variance 
effect 

L09 report a morphological merger fraction 

.C,Gs(^ = 0.6, Ms < -20) = OMStZn d D 

in the Groth field (open diamond in Fig.|7]i. How does this value 
compare with the one obtained in GOODS-S? If we use the same 
selection as in L09, this is, Mg < -20 galaxies with 0.35 < z < 
0.85, the major merger fraction in GOODS-S is 

.CcooDsfe = 0.6, Mb < -20) = 0.032!°;°'^. (12) 

We can see that both values are consistent within their errors. 
Because both values are determined using the same methodol- 
ogy and sample selection, the difference of A/^ = 0.013 may 
be explained by cosm ic variance, denoted by o-y. Following 
ISomerville et al.l (|2004|) . we infer that the effect of cosmic vari- 
ance for the typical merger density (~ 10 '* Mpc"^, Sect. 16.4b 
and GOODS-S/Groth volume is ~ 60%. That is, we expect 
o-v ~ 0.027 in Groth and ct^ ~ 0.019 in GOODS-S, so the dif- 
ference between both merger fraction determinations can indeed 
be explained as cosmic variance effect. Averaging both values, 
the morphological merger fraction at z = 0.6 is 

f7^(z = 0.6, Mb < -20) = 0.038 + 0.012, (13) 

where the error is t he expected cry ~ 30% due to combining two 
separate fields (see lSomerviUe et aLll2004 for details). 

6.2. Morpfiological merger fraction evolution in previous 
studies 

In Sect. 15.11 we obtained the values of m and /"''''(O) that de- 
scribe the morphological merger fraction evolution in GOODS- 
S. In this section we compare these values with those in the 
literature, where morphological works in Z?-band selected sam- 
ples are common. L09 study the merger fraction for Mb < -20 
galaxies in Groth by asymmetries and taking into account the 
experimental error bias. Combining their results with the litera- 
ture, they obtain m = 2.9 + 0.8, consistent to within ~ Icr with 
our result. iLotz et akl (^008a) study the merger fraction in an 
Mb < -18.83 - 1.3z selected sample by G and M20 morpho- 
logical indices. Their results alone suggest m = 0.23 + 1.03, 
but when combined with others in the literature they obtain 
m = 2.09 + 0.55. The first case does not match the local mor- 
phological merger fraction by De Propris et al. (2007): with a 
similar luminosity cut, Mb < -19, and taking into account the 
different methodologies (see L 09, for details), the mer ger frac- 
tions are very diff erent, 0.006 dDe Propris et al.ll2007h vs 0.07 
(iLotz et al. 2008a). Because of this, the second m value is pre- 
ferred. Kampczyk et al. (2007) study the fraction of visually dis- 
torted galaxies in SDSS0(S1 oan Digital Sky Survey, local value) 
and COSMOS (z ~ 0.7 value) for Mb < -19.15 galaxies. They 
find that m = 3.8 ± 1 .8, higher than our va lue, but consistent to 
within ~ Icr. Finallv. C onselice et al.l ("2003) study the morpho- 
logical merger fraction of Mb < -20 by asymmetries. However, 
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Fig. 7. Morphological merger fraction vs redshift for Mb < -20 (open squares) and M* > lO'" M© galaxies (bullets). The error 
bars do not include cosmic variance (Sec. 16.11 ). The grey/black solid lines are the least-squares fit of /"''^ (z) = /"''''(0)(1 + z)"' to 
the data in the luminous/mass case, respectively. The open tri angle is the Mb < -20 estimate by L09 of the lDe Propris et aTl (l2007h 
Mb < -19 merger fraction, open inverted triangles are fromjLotz et al.l (I2008al), open diamond is for Mb < -20 galaxies in Groth 
strip from L09, crosses are for M^, > 10'° Mq ga l axies from[Conselice et al. (2009), and filled pentagons are minor+m ajor mergers 
forM* ^ 5 X 10'° Mo galaxies from lJogee et al](l2009l) . The dashed lines marks the major merger fraction expected by! Jogee et al.l 
(120091) . 



due to the small area of their survey, they have high uncertain- 
ties in the merger fraction at z < 1, so we do not compare our 
results with theirs. In summary, the morphological major merger 
fraction evolution in Mb samples up to z ~ 1 is consistent with 
a OT = 2.2 + 0.4 evolution (weighted average of the previous m 
values), although more studies are needed to understand its de- 
pendence on different luminosity selections. 

The only previous morphologic al merger fractions in 
10'° M q selected samples are from lConselice et alj ( l2003ll2008l 
20091). The small a reas in the f irst two studie s m 
Conselice et 311120031 and UD10 in lConsehce et al.ll2008l) make 
their z < 1 values highly undetermined, and we use their z > 1 
values to constrain the merger fraction evolution at higher red- 
shifts in Sect. |631 IConseUce et al.1 (12009) find m = 3.8 ± 0.2. 
This value is lower than ours, but it is higher than typical values 
in B-band studies, supporting the hypothesis that merger frac- 
tion evolution in mass-selected samples is more important than 
in luminosity-selected samples. 

Other asym metry studies have us ed different selection cri- 
teria from ours: ICassata et al.l (l2005h obtain a merger fraction 
evolution m - 2.2 ± 0.3 in an niK^ < 20 selected sample, and 
combining their results with others in the literature. [Bridge et alj 



^ http://www.stsci.edu/ftp/science/hdf/hdf.html 
http://www.stsci.edu/hst/udf 



(l2007h perform their asymmetry study on a 24jum-selected sam- 
ple (L/R > 5.0 X 10'° Lq), finding m - 1.08. However, these 
values are difficult to compare with ours because studies with 
selection s in differ ent bands yields different results dBundv et al.l 
2004;,Rawatet al.|[2008; L09). 



6.3. Merger fraction evolution at iiigiier redsliift 

Merger fraction studies of Mj, > 10'° Mq g alaxies at redshift 
higher than z ~ 1 are rare. [Ryan et akl (l2008h address the p rob- 
lem with pair statistics, while Conselice et alj (I2003L 1200 8*) use 
asymmetries. Both these studies conclude that the merger frac- 
tion shows a maximum at z > 1 .5 and decreases at higher z. This 
tells us that we cannot extr apolate the power-law fi t (Eq. ifTOl ) to 
high redshift. Fortunately, IConselice et al.l (l2008l) perform their 
study by asymmetries, providing us with a suitably high redshift 
reference. Note that, although Conselice et al. (2008) treated the 
loss of information with redshift, they do not take into account 
the overestimat i on due to the experimental errors. Because the 
IConselice et al.l (l2008l) study is performed in UDF, which is 
located in the GOODS-S area, we apply a 0.5 factor to the 
(Conselice et al. (2008) merger fractions based on the results of 
Section|5] In Fig.[8]we show the corrected lConselice et al.l(l2008l) 
merger fractions (white dots) and our data (black dots). Note that 
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Fig. 8. Morphological merger fraction vs redshift fo r > 
10^° Mq galaxies. Data are from lConseUce et aTl (l2008l open cir- 
cles) and this work (bullets). The black solid line gives the least- 
squares power-law fit to our data, fm^'^iz) = 0.00 1(1 -Hz)^'^, while 
the dashed grey hne is the least-squares fit to IConselice et alj 
(I2008h data at z > 1 .2 and ours at z < 1 .2, f7^(z) = 0.00034(1 + 

^y0.5g-0.57(l+i)2 



the previous power-law fit to our data (black solid line. Sect. 15. U 
fails to explai n the merge r fraction values at z > 1.5. 

FoUowing lConselica (l2006b) . we parameterize the observed 
tendency as 



/,7\z) = a(l+z)'"e/^('"* 



(14) 



where the local merger fraction is given by /^''^(O) - aexpiP). 
This form is also obtained for the evo lution of the m erger rate on 
the basi s of Press-Sch echter theory ( Carlbergll990l) . The best fit 
of the Conselice et al.l ( 12008) data at z > 1.2 and ours at z < 1.2 
yields a = 3.4 X 10^^*, m = 10.5, and jS = -0.57 (gray dashed 
line in Fig. [8]l. With these value s the merger f raction peaks at 
Zpeak - 2, in good agreement with IConselice et al. (2008). 

The previous parameterization implies that the merger frac- 
tion drops at z > Zpe a k - 2 ; being ~ 0.01 at z ~ 4. On the other 
hand. lHopkins et aTl (l2008h models suggest that the merger frac- 
tion of Mi, > lO'" Mq galaxies still grows at z > Zpeak, being 
~ 0.30 at z ~ 4. In fact, the data in Fig.|8]can also been fitted by 



.mph. . _/0.001(l+z)^-4 Z<Zc 
( Jm.c 4. ^ 4.C 



(15) 



where /"c*" is a constant, and Zc is the redshift in which the 
merger fraction behaviour changes; that is, when 0.001(1 -H 
Zc)^'^ = /ill!?- With the two high redshift points in Fig. |8] we 
estimate that f^f -0.18 + 0.04, which yields Zc = 1 .62. Further 
studies are needed to constrain the merger fraction evolution at 
high redshift, but it is clear that the potential approximation is 
only valid at z 5 1.5. 

6.4. The major merger rate evolution 

We define the major merger rate (^RU)'''') as the comoving number 
of major mergers per Gyr within a given redshift interval and 
luminosity or stellar mass range: 



n7\z,M)^piz,M)f7\z,M)T-^]^, 



(16) 



Fig. 9. Morph ological merger rate vs redshift for Mb < -20 
(open triangle, iDe Propris et al.ll2007l: open squares, this work) 
and Mi, > 10'° Mq galaxies (bullets). The grey/black solid line 
is the least-squares fit ^R7\z) = iRr\0)(l + z)" to the data 
from the luminosity/mass-selected sample, respectively. 



where M = Mb I Mi, denotes the selection of the sample, 
p(z,M) is the comoving number density of galaxies at red- 
shift z brighter/more massive than Mb I Mi,, and is the 
merger timescale in Gyr for the asymmetry crite rion. To ob- 
tain p(z,Mb) we assume the Faber et al] (|2007|) luminosity 
function parameters, while to obtain n(z, M*) we assume the 
iPerez-Gonzalez et al.l (l2008l) mass function parameters. In ad- 
dition, we take T^ ^ = 0.35 - 0.6 Gyr. The lower value is from 
Conselice (2006b, N-body major merger simulations), and the 
higher from Lotz et a l. (2008b;, N-body/hydrodynamical equal- 
mass merger simulations). 

We summarize the merger rates in Table |5] and show these 
values in Fi^.0 white sy mbols a re for Mb < -20 galaxies (white 
triangle. iDe Propris et aLli2007t white squares, this work) and 
black dots for Mi, > lO'" Mq, while the grey/black solid line is 
the least-squares fit of ^R7\z) = !R7''(0)(1 +z)" function to the 
data in the luminosity-/mass-selected sample. The parameters of 
these fits are also summarized in Table |5] In spite of the very 
different merger fraction evolution, the merger rate evolution of 
both samples are similar: n = 3.3 ±0.8 in the luminosity sample, 
while « = 3.5 ± 0.4 in the mass sample. As in the merger frac- 
tion case, the results are not affected by the LSS (Sect. |572] i. The 
reason why the very different merger fraction evolution tuns into 
a similar merger rate evolution is the evolution over cosmic time 
of the number density of galaxies. The number of Mb < -20 
galaxies decreases by a factor 3 from z - 1 to z = 0, while the 
number of > 10'° Mq galaxies increases by a factor 3 in the 
same redshift range. 

We can compare our infer red merger rate w ith the post- 
starburst (PSB) rate reported bv Wild et af] (12009 *). The light of 
PSB galaxies is dominated by A/F stars. Such galaxies are iden- 
tifiable by their strong Balmer absorption lines compared to their 
mean stellar age as measured by their 4000 A break strength. 
PSB spectra indicate that the formation of O- and early B-type 
stars has s uddenly ceased i n the galaxy. The simulations per- 
formed by iJohansson et al.l (|2008|) find that the PSB phase can 
only be reached by disc-disc major merger remnants, so the 
PSB rate and our merger rate may be similar if an evolution- 
ary path connects both populations. The PSB rate, in the range 
0.5 < z < 1 and for M* > 10'° Mq galaxies (Salpeter IMF), 
is ^RpsB = (1.6-2.9) X 10""* Mpc"^ Gyr"', where the interval re- 
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Table 5. Major merger rates in GOODS-S 



Sample selection 


z = 0.4 " 


z = 0.725 " 


z = 0.975 " 


^K„,(0) 


n 


b 


Mb < -20 


1.2+'-;' 




-t q+l .') 


0.40 ±0.14'^ 


3.3 ± 


0.8 


> 10'" Mo 


9+^-^ 




9 Q+1.6 
"^■^-1.0 


0.29 + 0.06 


3.5 H 


I 0.4 



" In units of 10"" Mpc"' Gyr"'. 
Best ^'Tiz) = 'K 7\0)il + zY fit to the data. 
This fit includes the lDe Propris et aTl l l2007b local value. 



fleets the u ncertainty in the PSB phase time-scale (0.35-0.6 Gyr, 
IWild et all 12009). This value compares well with the inferred 
disc-disc major merger rate at that range, - (1.2-3.0) 

X 10 Mpc"^ Gyr Although the uncertainties in both studies 
are important, the result suggests that SPB galaxies can be the 
descendants of our distorted, disc-disc major merger remnants. 

6.4.1 . Number density of merger remnants 

If we integrate the merger rate over cosmic time, we obtain the 
number density of galaxies that have undergone a disc-disc ma- 
jor merger (prem) in a given redshift range: 



PKm(Zl,Z2) 



r 



n7\o)(i + zr 



dz 



HoEiz) 



(17) 



where E(z) = -\/^a + ^m(^ + z)'^ in a flat universe. We make 
this study only for the mass-selected sample because we can as- 
sume that stellar mass is additive: M^^izi) > M^,(z2) always for 
zi < Z2, and prein(0,z) is representative of the number density 
of local galaxies that have undergone a disc-disc merger since 
redshift z- The same cannot be said for the luminosity-selected 
sample: here the number density of objects above a given ab- 
solute magnitude can decrease with time, as it is not generally 
flie case that Mb(zi) < Msizi) for zi < zi- Using Eq. ( fTTI l for 
Prem(0, z) would Overestimate the number of local galaxies that 
have undergone a merger 

Comparing prem with the number of > 10'" Mq galaxies 
at redshift z\ , p(zi), we obtain the fraction of merger remnants. 



/rem(Zl,Z2) 



Prem(Zl,Z2) 

P(zi) 



(18) 



Applying Eq. ( fTSl l with the merger rate parameters of the 
mass sample from Table |5] and the mass functions from 
iPerez-Gonzalez et all (12008!). we obtain /rem(0, 1) = 8!^%. This 
is a low value that increases to /rem(0, 1.5) - 15;';g%. We 
take Z2 = 1-5 as an upper limit because our merger frac- 
tion parameterization is vah d to this redshift (Sect. I6.3l l and 
iPerez-Gonzalez et alj (l2008h mass functions are complete for 
> 10'° Mq galaxies also up to z ~ 1.5. Interestingly, we 
infer that /|-em(l-0, 1.5) - 21^]^^%, which is compatible with 
the fraction of bulge-dominated galaxies (E/SO/S a) at z ~ 1 
(Lope z-Sanjuan et al., in prep). The pair study of lBundv et alj 
j2009l) in GOODS-S reports /rem in the range 0.4 < z < 1.4. 
For M* > 2 X 10'" Mq galaxies they estimate /rem(0.4, 1.4) - 
15%-18%, which is in good agree ment with our inferre d value, 
/rem(0.4, 1.4) ~ 17%. Givcn that the lBundv et"ai](l2009l) studv is 
also sensitive to mergers between spheroids, and the mass limit 
in both studies is different, the quantitative agreement is remark- 
able. 

The most important error source in our results is the uncer- 
tainty in the lower redshift bin, especially in the mass-selected 



sample. We repeat our study with a higher merger fraction in 
this bin by a factor of two, /^''^(z — 0.4) = 0.012, and three, 
fm^^iz = 0.4) = 0.018. With these assumptions /rem(0, 1) in- 
creases to 12% and 18%, respectively. These values remain low, 
so our conclusions do not change. 

6.4.2. Number of mergers per massive galaxy 

As a complement to the previous section we calculate the num- 
ber of expected disc-disc major mergers per > lO'" Mq 
galaxy in a given redshift range. 



NMi.zi) 



r 



'^'n7\z) dz 



P(z) 



f7\o)(i+zr-' 



HoE(z) 
dz 



T^^aHqE{z) 



(19) 



Taking Eq. ( fT4b as the merger fraction parameterization we ob- 
tain A^m(0,3) = l-2!o2- In addition, we also obtain A^m(l,3) = 
1.0^Q2, with only 0.2 disc-disc major mergers in < z < 1. 
The results are the same if we take Eq. ( fTSl ) as the merger 
fraction parameterization. T he previous ff m(0, 3) value is lower 
than that inferred by Bluc k et al.l (l2009ln for > 10" Mq 
galaxies, A^in(0, 3) - l-8^y 4, which implies that more massive 
galaxies have higher number of mergers than less massive ones. 
On the other hand , our value is ~ 4 times lower than that of 
IConsehcd ( l2006bl) : IConselice et al] ( l2008i) . A^gm(0, 3) ~ 4.4 for 
Mi, > 10'" Mq galaxies. We suspect that two factors con- 
tribute to their high value. First, they use the galaxy merger frac- 
tion (/gm, the fraction of galaxies undergoing mergers) to obtain 
Ngmizx^zi), the mean number of galaxies that merge since zi to 
obtain a zi galaxy. This is roughly 2 times higher than (the 
number of merger events per galaxy). And second, they use clas- 
sical counting statistics which, as shown in Sectsl5]and l6.3l leads 
to an overestimate of the merger fraction by another factor of 2. 

These results suggest that most of the disc-disc merger ac- 
tivity of M* > 10'" Mq galaxies happened before z ~ 1, 
this kind of merger being important in galaxy evolution down 
to this redshift. It is important to recall that our methodology 
cannot detect spheroidal major mergers, so the role of these 
mergers in the evolution of the red sequence since z ~ 1 
(lBelIetal.ll2004tlFaber et al.ll2007l) cannot be addressed by our 
study. However, due to the paucity of spheroid al systems at 
z > 1 .2 (Cons elice et al.ll2005 ; ICassata et al.ll2005l) . one expects 
spheroidal major mergers t o be important at lower redshifts; i.e. 
z < 1.2. The simulations of lKhochfar & Silkl (|2008|) are m agree- 
ment with this picture: they find that the dry merger rate is two 
orders of magnitude less than the wet merger rate at z ~ 1.5, 
while they are similar at z ~ 0.6. In addition, they find that 



' We apply Eq. M9\ to their best power-law fit of the merger fraction 
and assume, as previously, that r,^,A = 0.35 - 0.6 Gyr. 
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Fig. 10. Comparison between the observed and simulated merger 
fractions for galaxies with > 10'° Mr. The observations are 
from this work (black bullets), and from Conselice et aTl (l2009l 
crosses). The predictions (gray solid lines ) are for major merg- 
ers w ith a time-scale of T^a = 0.4 Gyr (iBertone & Conselice! 
l2009h . 



the wet merger rate has its maximum at z ~ 1.3, and then de- 
clines by an order of magnitude until z = 0, a similar evolution 
to our results. To check these ideas 

we explore the relative importance of disc-disc mergers in the 
structural evolution of > 10'*' Mq galaxies in a forthcoming 
paper. 



6.5. Comparison with model predictions 

The comparison of the predicions by cosmological simulations 
with our results is not straighforward because we only detect 
disk-disk (i.e. wet) major m ergers, an d we select by stellar mass: 
the simulations from Stewart et al.l (12008) point out that the 
merger fraction depends on merger definition (minor vs major), 
selection criteria (halo mass, stellar mass, or luminosity) or the 

assumed merger time-scale^ 

The study from B ertone & Conseli^ (1200 9) provides pre- 
dictions for major mergers of M*. > lO'" Mq galaxies, assuming 
a merger time-scale of T^ a - 0.4 Gyr. In figurelTOlwe show the 
predictions (gray solid line), and the observational data from this 
work (black bullets) and Conselice et al. (2009, crosses). The 
pre dictions are in good ag reement with our observations, while 
the IConsehce et al.l (l2009l) values are higher than predicted by 
factors 2 to 6. However, this agreement must be taken as qual- 
itative more than quantitative because (i) cosmological simula- 
tions might underestimate th e major merger fract i on at that stel- 
lar mass, as pointed out by Bertone & Conseli^ (l2009l) . (ii) the 
predictions are for total (i.e. wet + dry) major mergers, while 
we report wet major mergers. This can lead i n a ~1% increase 
in the m erger fractions due to dry mergers ( Belletal.ll2006i 
iLotz et al. 2008b). And, (iii) the GOODS-S merger fractions 
might be lower that the cosmological value due to cosmic vari- 
ance (Sect. |6T| i. Despite these caveats, the agreement is remark- 
able. 

On the other hand, the simulations of IStewart et al.l (l2009l) 
provide /rem(0, 2) for major wet mergers in M* > lO'" Mq galax- 
ies: they predict /rem(0, 2) ~ 10%-20%, in good agreement with 
our inferred /,-eni(0, 2) -21+''%. Finallv. IWeinzirl et al.1 (l2009t) 
compare their study of M^, > 10'" Mq local s pira l galaxies with 
the pred ictions by the lKhochfar & SiUd (l2006l) and lHopkins et all 
(l2009bl) models. They find that both models are able to explain 



the observed bulge-to-total ratio (BIT) distribution, and predict 
that only 13-16% of today's BjT < 0.75 spirals have undergone 
a major merger since z = 2. If we assume that all the disc-disc 
major mergers since z = 2 have enough g as to re-form a disc in 
the merger remnant dHopkins et alj20()9al) . our/rem(0, 2) ~ 20% 
value is an upper limit to the models' predictions, so both are 
compatible with our results. 

7. CONCLUSIONS 

We have computed the disc-disc major merger fraction and its 
evolution up to z ~ 1 in the GOODS-S field using morphologi- 
cal criteria. We quantify and correct for the bias due to varying 
spatial resolution and image depth with redshift by artificially 
redshifting the galaxy images to a common reference redshift of 
Zd = 1 • More importantly, we successfully account for the spill- 
over of sources into neighbouring bins caused by the errors in 
asymmetry indices and in Zphot, through the use of an ML method 
developed by LGB08. In every case we obtain merger fractions 
lower than 0.06, in agreement with the merger fraction deter- 
mination for the Groth field (L09). The main improvement in 
our study over previous determinations is the robust methodol- 
ogy that takes into account the signal-to-noise variation of galax- 
ies with z and the observational errors: previous morphological 
studies using classical counts overestimate the disc-disc major 
merger fractions by factors of ~ 2. 

The merger fraction evolution in luminosity- and mass- 
selected samples are, respectively. 



/r'fcMB< -20) = 0.013(1 +z) 



1.8 



/r\z,M* > 10'" Mo) = 0.001(1 + z?-\ 



(20) 
(21) 



We study the effect of the LSS on these results and find that 
merger fractions do not change substantially. 

When we compute the merger rate for both samples, the 
very different merger fraction evolution becomes a quite simi- 
lar merger rate evolution: 



^R7\z) = 0.40 X 10"'*(1 + zt"" Mpc"-*Gyr 
for Mb < -20 galaxies and 



^R7\z) = 0.29 X 10"'*(1 + zy-" Mpc^-'Gyr 



(22) 



(23) 



for M* > 10'" Mq galaxies. This similar evolution is due to the 
different number density evolution with redshift: the number of 
> 10'" Mq galaxies increases with cosmic time, while the 
number of Mg < -20 galaxies decreases. 

The previous merger rates imply that only ~ 8% of today's 
Mi, > 10'" Mq galaxies have undergone a disc-disc major 
merger since z ~ 1. Interestingly, ~ 21% of these galaxies at 
z ~ 1 have undergone a disc-disc major merger since z ~ 1.5, 
which is compatible with the fraction of bulge-dominated galax- 
ies (E/SO/Sa) at z ~ 1 (Lopez-Sanjuan et al., in prep). This sug- 
gests that disc-disc major mergers are not the dominant process 
in evolution of > 10'" Mq galaxies since z ~ 1, with only 
0.2 disc-disc major mergers per galaxy, but may be an important 
process at z > 1, with ~ 1 merger per galaxy at 1 < z < 3. 

The most important error source in these results is the uncer- 
tainty in the lower redshift bin, especially in the mass-selected 
sample. More studies are needed to improve the statistics at low 
redshift and avoid cosmic variance effects. Another important is- 
sue is the sample definition, given that merger fraction depends 
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on mass and luminosity: larger samples permit us different selec- 
tion cuts in luminosity and mass, thus improving our knowledge 
of the importance of disc-disc major mergers in galaxy evolu- 
tion. 
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